(§)crossMark rescarch papers 

BiobgS^^^ structures of cytosolic and plastid-located 

Crystallography glutamine Synthetases from Medicago truncatula 



ISSN 1399-0047 



reveal a common and dynamic architecture 



Eva Torreira,''* Ana Rita 
Seabra,''* Hazel Marriott,'^ Min 
Zhou,'^ Oscar Llorca,* Carol V. 
Robinson,'^ Helena G. Carvalho,'' 
Carlos Fernandez-Tornero^* and 
Pedro Jose Barbosa Pereira''* 



''Chemical and Physical Biology, Centre de 
Investigaciones Biologicas - CSIC, Ramiro de 
Maeztu 9, 28040 Madrid, Spain, ''Molecular 
Biology of Nitrogen Assimilation, IBMC - 
Institute de Biologia Molecular e Celular, 
Universidade do Porto, Rua do Campo Alegre 
823, 4150-180 Porto, Portugal, ""Chemistry 
Research Laboratory, University of Oxford, 
South Parks Road, Oxford OXl 3QZ, England, 
and ''Biomolecular Structure Group, IBMC ~ 
Institute de Biologia Molecular e Celular, 
Universidade do Porto, Rua do Campo Alegre 
823, 4150-180 Porto, Portugal 



t These authors share first authorship. 



Correspondence e-mail: cftornero@cib.csic.es, 
ppereira@ibmc.up.pt 



The first step of nitrogen assimilation in higlier plants, the 
energy-driven incorporation of ammonia into glutamate, is 
catalyzed by glutamine synthetase. This central process 
yields the readily metabolizable glutamine, which in turn is 
at the basis of all subsequent biosynthesis of nitrogenous 
compounds. The essential role performed by glutamine 
synthetase makes it a prime target for herbicidal compounds, 
but also a suitable intervention point for the improvement 
of crop yields. Although the majority of crop plants are 
dicotyledonous, little is known about the structural organiza- 
tion of glutamine synthetase in these organisms and about the 
functional differences between the different isoforms. Here, 
the structural characterization of two glutamine synthetase 
isoforms from the model legume Medicago truncatula is 
reported: the crystallographic structure of cytoplasmic GSII- 
la and an electron cryomicroscopy reconstruction of plastid- 
located GSII-2a. Together, these structural models unveil a 
decameric organization of dicotyledonous glutamine synthe- 
tase, with two pentameric rings weakly connected by inter-ring 
loops. Moreover, rearrangement of these dynamic loops 
changes the relative orientation of the rings, suggesting a 
zipper-like mechanism for their assembly into a decameric 
enzyme. Finally, the atomic structure of M. truncatula GSII-la 
provides important insights into the structural determinants 
of herbicide resistance in this family of enzymes, opening new 
avenues for the development of herbicide-resistant plants. 
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1. Introduction 

Glutamine synthetase (GS; EC 6.3.1.2) is a vital enzyme in the 
nitrogen metabolism of all living systems. GS catalyzes the first 
step in nitrogen assimilation, the ATP-dependent ligation of 
ammonia and glutamate to yield glutamine (1), which is then 
used for the biosynthesis of essentially all nitrogenous 
compounds in the cell (Lea & Miflin, 2010), 

GS 

glutamate -|- ammonia -I- ATP 

glutamine + ADP -h Pj -h HjO. (1) 

The assimilation of ammonia is important for all organisms, 
but it is of special interest in plants, both because it is the main 
entry route of inorganic nitrogen into organic composition and 
because nitrogen is a major limiting factor in plant growth and 
productivity. The reaction catalyzed by GS is a key component 
of plant nitrogen-use efficiency (NUE) and plant yield, and 
an increase in GS activity can potentially enhance plant 
productivity (Lea & Miflin, 2010). Conversely, GS inhibition 
leads to plant death, and a number of herbicides have been 
designed on this principle. An in-depth understanding of the 
molecular details of glutamine synthetase activity is therefore 
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of crucial importance not only to improve NUE but also for 
the design of herbicide-resistant plants. 

Glutamine synthetase is a fascinating, complex and highly 
regulated enzyme found in all living organisms (Kumada et al. , 
1993), but GS proteins differ significantly in molecular mass, 
sequence and quaternary structure among the different forms 
of life. Three major classes, referred to as GSI, GSII and 
GSIII, have been described (Woods & Reid, 1993). GSI, the 
best-characterized class, includes enzymes typically found in 
prokaryotes that display a homododecameric arrangement of 
subunits with masses ranging from 44 to 60 kDa. Following 
the groundbreaking determination of the atomic structure of 
GlnA from Salmonella typhimurium by X-ray crystallography 
(Almassy et al, 1986), the structures of GSI from Myco- 
bacterium tuberculosis (Gill et al, 2002) and Synechocystis sp. 
(PDB entry 3ng0; L. Saelices, D. Cascio, F. J. Florencio & M. I. 
Muro-Pastor, unpublished work) have greatly contributed to 
current knowledge of the structural features and enzymology 
of GS. Typical eukaryotic GSs belong to class II, forming 
homodecamers of 35-50 kDa subunits. The first three- 
dimensional GSII structure determined was that from maize 
(Unno et al, 2006), followed by the mammalian (Krajewski et 
al, 2008) and yeast (He et al, 2009) enzymes, all of which are 
present in the cytoplasm. Finally, GSIII has been identified in 
cyanobacteria (Reyes & Florencio, 1994) and two unrelated 
anaerobic bacteria (Southern et al, 1986; Goodman & Woods, 
1993). Although its 75-83 kDa subunits were thought to be 
organized as homohexamers (Reyes & Florencio, 1994), the 
crystal structure of Bacteroides fragilis GSIII revealed a 
dodecameric architecture showing significant structural 
differences from GSI (van Rooyen et al, 2011). Interestingly, 
genes encoding prokaryotic-type GSI-like molecules have also 
been found in plants (Mathis et al, 2000; Doskocilova et al, 
2011) and mammals, and display a structural organization 
similar to their bacterial homologues, as confirmed by the 
electron-microscopy structure of mouse lengsin ( Wyatt et al , 
2006). On the other hand, a few soil-dwelling bacteria 
(belonging to the Rhizobiaceae, Frankiaceae and Strepto- 
mycetaceae families) contain an eukaryotic-type GS (Carlson 
& Chelm, 1986). Whereas the function of these proteins is still 
unclear, they probably result from lateral gene transfer during 
evolution (Ghoshroy et al, 2010). 

The overall geometry of the active site is the most 
conserved structural element amongst GS enzymes (Eisen- 
berg et al, 2000; Unno et al, 2006; van Rooyen et al, 2011). 
This funnel-shaped pocket is formed by the N-terminal 
j6-grasp domain of one monomer and the highly curved 
j6-sheet in the C-terminal catalytic domain of the neighbouring 
monomer. The pocket closes upon nucleotide binding and 
several additional changes occur upon glutamate attachment 
(Krajewski et al, 2008). Despite conservation of the active-site 
architecture, eukaryotic GSII does not contain the adenyla- 
tion loop which regulates GS activity in bacteria. Another 
remarkable structural difference between GSI and GSII is the 
inter-ring association. While hexameric GSI rings tightly 
interact by insertion of the C-terminal helix of each monomer 
into a hydrophobic cavity on the corresponding subimit of the 



opposite hexameric ring, GSII rings interact weakly and only 
through internal loops. 

In plants, the ammonia for GS activity may be derived from 
primary sources (nitrogen fixation, nitrate assimilation and 
ammonia uptake) or from secondary metabolism such as 
photorespiration and amino-acid catabolism; therefore, the 
enzyme plays essential roles in both primary nitrogen assim- 
ilation and nitrogen recycling (Bernard & Habash, 2009; Lea 
& Miflin, 2010). Consistent with the diversity of metabolic 
roles, GS exists in plants as a number of isoenzymes, located 
both in the cytosol (GSII-1) and in the plastids (GSII-2), that 
are encoded by a small gene family. Each of the GS genes 
displays distinct patterns of expression in different organs of 
the plant and their encoded enzymes participate in different 
metabohc processes. In the model legume Medicago trunca- 
tula, the GS gene family consists of four expressed genes, 
MtGSII-la, MtGSII-lb, MtGSII-2a and MtGSH-lb (Stanford 
et al, 1993; Carvalho, Lescure et al, 2000; Carvalho, Lima et 
al, 2000; Melo et al, 2003; Seabra et al, 2010), the latter of 
which is unique to M. truncatula and closely related species. 
The other three GS genes are expressed in almost all organs of 
the plant, but in a cell-specific manner. MtGSII-la is expressed 
in the vascular bundles of all plant organs (Carvalho, Lescure 
et al, 2000; Carvalho, Lima et al, 2000) and is highly up- 
regulated in the nodules, accounting for the production of over 
90% of the total nodule GS activity. The M^G5/7-7a-encoded 
isoenzyme is responsible for the assimilation of ammonia 
released by nitrogen fixation (Carvalho, Lima et al, 2000). The 
enzyme encoded by MtGSII-lb appears to play a house- 
keeping role, being expressed at lower levels and almost 
everywhere in the plant, whereas MtGSII-2a is highly 
expressed in all photosynthetic tissues and is mainly respon- 
sible for the reassimilation of ammonia released by photo- 
respiration (Carvalho & Cullimore, 2003). Finally, MtGSII-2b 
is exclusively expressed during seed filling and its physiolo- 
gical role is probably therefore related to protein storage 
(Seabra et al, 2010). 

The specialization of closely related isoenzymes to a parti- 
cular metabolic context certainly reflects specific functional 
determinants that make them adapt better to each environ- 
mental condition. However, very little is known regarding the 
structural features and catalytic determinants of the different 
isoenzymes. The only plant GS enzyme structurally char- 
acterized to date is maize cytosolic ZmGSII-la, which proved 
to be a decamer (Unno et al, 2006). However, an octameric 
architecture has been suggested for cytosoHc GSII from 
Phaseolus vulgaris (Llorca et al, 2006; Betti et al, 2012), while 
the structure of plastid-located GS is unknown. Comparison of 
experimental models of closely related proteins can assist in 
the identification of critical structural components responsible 
for enhancement of activity and in the design of herbicide- 
resistant enzymes. 

Here, we present the unique crystallographic structure of a 
cytosolic GS from a dicotyledonous plant, the model legume 
M. truncatula. We also present the electron cryomicroscopy 
(cryo-EM) structure of its plastid-located GS enzyme and 
demonstrate that it similarly organizes in homodecamers 
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Table 1 

Data-collection and refinement statistics. 

Values in parentheses are for the outermost resolution shell. 

Data collection 



Space group 


P2i 


Unit-cell parameters (A, 


a = 99.3, b = 101.7, c = 188.1, 




^ = 103.7 


Resolution range (A) 


101-2.35 (2.48-2.35) 


Reflections (measured/unique) 


935718/149655 


Completeness (%) 


99.0 (96.8) 


Multiplicity 


6.3 (5.6) 




0.077 (0.338) 


Mean {//<t(/)> 


19.2 (4.0) 


Monomers per asymmetric unit 


10 


Matthews coefficient (A' Da~^) 


2.2 


Solvent content (%) 


44.7 


Refinement 




Resolution range (A) 


96.5-2.35 


R factorf/free R factor§ (%) 


17.5/21.7 


Unique reflections (working/test set) 


141942/7518 


Water molecules 


969 


Total No. of atoms 


24909 


R.m.s.d., bond lengths (A) 


0.008 


R.m.s.d., bond angles (°) 


1.165 


Ramachandran plot statistics 




Residues in allowed regions (%) 


100 


DPIH (A) 


0.21 



t J?merge = UmllMhkl) - (I(hkl))\IY.HuT.Mhkl), where I,(hkl) is the observed 
intensity and {IQikl)) is the average intensity of multiple observations of symmetry- 
related reflections. % R factor = Y,hu | l^'obs I - I ^'cak 1 1 /Y,hti I ^obs I . where |F„bsl and \F^\ 
are the observed and calculated structure-factor amplitudes, respectively. § The free R 
factor is the cross-validation R factor computed for a randomly chosen subset of 5% of 
the total number of reflections which were not used during refinement. T| Diffraction- 
data precision indicator. 

composed of two pentameric rings. We show that these rings 
associate weakly and propose a dynamic zipper-Uke 
mechanism for inter-ring assembly. A comparison between 
these structures leads us to identify common features in 
enzymes that function under different biochemical contexts 
and differ in catalytically properties, and provides insights into 
the structural determinants of herbicide resistance. 



2. Methods 

2.1 . Recombinant protein expression and purification 

The cloning procedures for N-terminally His-tagged 
proteins have been described previously (Seabra et al, 2009). 
The cDNAs coding for M. tmncatula MtGSII-la and MtGSII- 
Ib and the sequences coding for the mature MtGSII-2a and 
MtGSII-2b polypeptides were amplified by PGR from the 
vector pTrc99A (Carvalho et al, 1997; Melo et al, 2003; Seabra 
et al, 2010) using the forward primer 5'-TTGACAATTAAT- 
CATCCGGC-3' and the specific reverse primers 5'-TGG- 
TTGTGGTCGACTGGTTTCC-3' (MtGSII-la), 5'-AGCGT- 
GGTGTCGACTGGTTTCC-3' (MtGSII-lb), 5'-AATAGA- 
TGTCGACTTTCAATGC-3' (MtGSII-2a) and 5'-AATAGA- 
TGTCGACCTTCAATGC-3' (MtGSII-2b). The same forward 
primer, which annealed in the vector backbone, was used to 
amplify the four cDNAs. The reverse primers were designed to 
specifically abohsh the stop codon of each GS cDNA and to 
include the SaK restriction sequence. The PGR fragments 
were partially digested with the restriction enzyme Ncol and 



fully digested with SalV. The PGR digests were cloned into the 
Nco\ and Xho\ sites of pET-24-d-T vector, a derivative of 
pET-24-d(-H) (Novagen) containing a G-terminal His tag and 
a thrombin cleavage site. The resulting plasmids encode 
C-terminally Hisg-tagged GS fusion proteins in which the 
sequence LVPRGSVEHHHHHH foUows the GS coding 
sequences. The four expression constructs were sequenced. 

To express the His-tagged GS proteins, Escherichia coli 
BL21 CodonPlus (DE3) RIL cells (Stratagene) harbouring 
the pET-24d-GS plasmids were first cultured in LB medium at 
37°C to mid-exponential growth (ODgoo = 0-5). Expression of 
the recombinant proteins was then induced with IPTG (final 
concentration 1 mM) and cell growth continued overnight 
at 20° C. The cells were harvested by centrifugation at 2800g, 
resuspended in extraction buffer (10 mM HEPES buffer 
pH 7.4, 10 mM magnesium sulfate, 5 mM glutamate, 500 mM 
NaCl, 20 mM imidazole), disrupted by sonication and centri- 
fuged (60 min, 38 OOOg, 4°C) to remove cell debris. The crude 
protein extract was filtered through a 5 \im low-protein- 
binding filter and loaded onto a 5 ml Ni-Sepharose column 
(GE Healthcare) equilibrated with buffer A (10 mM HEPES 
buffer pH 7.4, 500 mM NaCl, 20 mM imidazole). Elution of 
the bound fusion protein was achieved by increasing the 
imidazole concentration in buffer A to 230 mM. The GS- 
containing fractions were pooled, dialyzed against 10 mM 
HEPES PH 7.4 and concentrated to 5 mg ml"^ on a centri- 
fugal concentration device with a 10 kDa molecular-weight 
cutoff membrane (Vivaproducts). 

2.2. Crystallographic structure determination and refinement 

M. tmncatula GSII-la was crystallized in space group P2i 
and a complete 2.35 A resolution X-ray diffraction data set 
was collected from a single cryocooled crystal on ESRF 
beamline ID14-EH2 as described previously (Seabra et al, 
2009). The data were processed with MOSFLM (Leshe & 
Powell, 2007) and scaled with SCALA (Evans, 2006) from the 
CCPA suite (Winn et al, 2011). The structure was solved by 
molecular replacement with Phaser (McCoy et al, 2007) using 
the coordinates of Zea mays GSla (Unno et al, 2006) as search 
model with all non-identical nonglycine residues truncated to 
alanine (Seabra et al, 2009). 

The model was iteratively improved by alternating cycles of 
manual model building with Coot (Emsley et al, 2010) and 
refinement with CNS (Brunger, 2007). The resulting isotropic 
model was further refined using TLS parameters and torsion- 
based NCS restraints as implemented in PHENIX (Afonine et 
al, 2012). The refinement statistics are summarized in Table 1. 
The final model has good stereochemistry, with 100% of the 
residues in the allowed regions of the Ramachandran plot. It 
comprises residues Metl-Tyr284 and Asn302-Pro356 in chain 

A, Metl-His277, Val315-Val317 and Tyr328-Leu353 in chain 

B, Metl-Lys278, Ser314-Val317 and Tyr328-Pro356 in chain 

C, Metl-Arg276, Ser314-Val317, Gly327-Asp331 and 
Met338-Pro356 in chain D, Metl-Glu259, Ile265-Arg276, 
Phe329-Asp331 and Met338-Pro356 in chain E, Metl-His277, 
Ser314-Val317, Gly327-Arg332 and Asn337-Pro356 in chain 
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Figure 1 

(a) Amino-acid sequence alignment of plant glutamine synthetases, including M. truncatula GSII-la (MtGSIIIa; UniProt 004998), GSII-lb (MtGSIIlb; 
UniProt 004999), GSII-2a (MtGSII2a; UniProt 084UC1), GSII-2b (MtGSII2b; UniProt E1ANG4) and Z. mays GSII-la (ZmGSIIIa; PDB entry 2d3a; 
Unno et al, 2006; differing from UniProt B9TSW5 by a I353V replacement), performed with ClustalW (Larkin et al, 2007). The transit peptides of 
M. truncatula GSir2a and GSII-2b were omitted. Red, green and white backgrounds are used for high, medium and low conservation, respectively. 
Numbers above the alignment and secondary-structure elements (a-helices in red, 3]o-helices in yellow and /6-strands in blue) correspond to the crystal 
structure of MtGSII-la described here. A dashed line denotes disordered residues in the crystallographic model. The ligand-contacting residues in 
Z. mays GSII-la are indicated below the alignment with red (ATP), blue (glutamate), green (ammonia) or open (metal coordination) triangles. A pink 
line above the alignment indicates the structurally distinct regions between MtGSII-la and ZmGSII-la. Prepared with Aline (Bond & Schiittelkopf, 
2009). (b) Inhibition of M. truncatula GSII by increasing concentrations of PPT (left) and MSO (right). The GS activity was normalized to that found in 
the absence of inhibitor. PPT and MSO were titrated using standard synthetase reactions with 100 or 20 mM glutamate, respectively. 



984 Torreira et a/. • Cytosolic and plastid-located glutamine synthetases 



Acta Cryst. (2014). D70, 981-993 



research papers 



F, Metl-His277, Ser314-Val317 and Gly327-Pro356 in chain 

G, Metl-His277, Tyr328-Arg332 and Asn337-Pro356 in cliain 

H, Metl-Arg276, Ser314-Val317, Tyr328-Asp331 and 
Met338-Pro356 in chain / and Metl-Arg276, Gly327-Arg332 
and Asn337-Pro356 in chain / of the glutamine synthetase 
decamer present in the asymmetric unit. There was no inter- 
pretable electron density for the remaining amino-acid resi- 
dues, including the totality of the N-terminal affinity tags, 
which were therefore excluded from the model. The refined 
coordinates and structure factors were deposited in the PDB 
(http://www.pdb.org) under accession code 4is4. 

2.3. Electron microscopy 

Purified GSII-2a (O.Smgml"^) was vitrified using glow- 
discharged QUANTIFOIL R 2/2 holey grids and a Vitrobot 
(Gatan) according to previously described methods (Dubo- 
chet et al, 1988). Frozen grids were transferred to a JEM- 
2200FS microscope at low temperature (— 170°C) using a cryo- 
transfer holder (Gatan). The microscope was operated under 
cryogenic conditions at an accelerating voltage of 200 kV. 
Data were collected imder low-dose conditions on a 4k x 4k 
CCD (Gatan) at a magnification of 83 393, with a pixel size of 

I, 62 A at the specimen level. The contrast transfer function of 
the micrographs was determined with ctfftnd3 (Mindell & 
Grigorieff, 2003). 7111 particle images were selected manually 
with boxer from the EMAN package (Ludtke et al, 1999). 
Reference-free two-dimensional alignment of the data set was 
performed with ML2D and CL2D from the Xmipp package 
(Scheres et al., 2005). Several starting models were generated 
from top-view and side-view averages applying D5 symmetry, 
as suggested by mass spectrometry, and were refined using 
the projection-matching algorithm in Xmipp. For the final 
reconstruction, only the 3117 particles corresponding to side 
views with parallel rings (see §3) were employed. The hand 
of the three-dimensional reconstruction was decided by 
comparison with available crystal structures of homologous 
proteins. The final resolution was estimated to be 20 A using 
Fourier shell correlation (FSC; 0.5 correlation-coefficient 
cutoff). Fitting of atomic structures into the cryo-EM density 
was carried out using UCSF Chimera (Pettersen et al, 2004). 

2.4. Native mass spectrometry 

A 20 |il aliquot of the MtGSII-2a complex was buffer- 
exchanged into ammonium acetate buffer pH 7.5 using Micro 
Bio-Spin 6 columns (Bio-Rad). To acquire MS spectra under 
nondenaturing conditions, 2 ^il ahquots were electrosprayed 
from gold-coated borosilicate capillaries prepared in-house. 
Spectra were recorded on a Q-Tof 2-type mass spectrometer 
(Waters) modified for high mass detection and adjusted to 
preserve noncovalent interactions (Hernandez & Robinson, 
2007). MS experiments were performed at a capillary voltage 
of 1700 V and a cone potential of 50 V. In experiments where 
high activation was employed, the collision energy was 
increased from 100 to 200 V. LC-MS analysis of the MtGSII-2a 
monomer was carried out on an NCS-3500RS nano-LC system 
(Thermo) equipped with a nano-UV detector set at 214 and 



280 nm. The MtGSII-2a sample was prepared in solvent A 
(0.05% trifluoroacetic acid; TFA) and 1 |il was applied onto 
a nano PS-DVB reverse-phase monolithic column (100 |xm 
internal diameter x 25 cm; Thermo) equilibrated with 90% 
solvent A and 10% solvent B (0.04% TFA, 90% acetonitrile). 
MtGSII-2a monomer was eluted with a linear gradient of 10- 
70% solvent B in 15 min at a flow rate of 600 nl min^'. The 
column effluent was delivered to a nanoflow ESI source and 
analyzed by MS on a QSTAR XL mass spectrometer (Applied 
Biosystems). 

2.5. Analysis of enzyme-activity inhibition 

The activities of the His-tagged GS isoenzymes (GSII-la, 
GSII-lb, GSII-2a and GSII-2b) were determined by quanti- 
fying the y-glutamyl hydroxamate produced by the synthetase 
reaction (GSs) as described previously (Culhmore & Sims, 
1981; CuUimore et al, 1982). The reactions started by the 
addition of 10 |ig purified His-tagged GS isoenzyme to the 
assay mixture (100 mM Tris buffer pH 7.8, 100 mM glutamate, 
8 mM ATP, 8 mM hydroxylamine, 16 mM MgS04) were 
allowed to proceed at 30°C for 15 min and were stopped by 
adding one volume of stop solution (0.37 M FeCls, 0.67 M 
HCl, 0.2 M TCA). GS specific activity was determined as the 
amount of /-glutamyl hydroxamate produced per minute per 
milligram of purified protein ((xmol min^^ mg~^). GS inhibi- 
tion was evaluated by adding increasing concentrations of 
either phosphinothricin (PPT; Fluka; 0, 0.005, 0.01, 0.05, 0.1, 
0.5, 1, 1.5, 2.5, 5, 10 and 20 mM) or methionine-j'-sulfoximine 
(MSG; Sigma; 0, 0.5, 1, 1.5, 2.5, 5, 10, 20, 40 and 60 mM) to the 
reaction mixture. GS inhibition assays were performed at two 
different glutamate concentrations (20 and 100 mM). IC50 was 
calculated by nonlinear regression using the LogIC50 equa- 
tion of GraphPad PrismS (GraphPad Software). 

3. Results 

3.1. Differential inhibition of M. truncatula GSII isoenzymes 

The four M. truncatula isoenzymes are highly conserved, 
displaying 70-94% sequence identity (Fig. la). Two small 
extensions (seven and 16 residues long at the N-terminus and 
the C-terminus, respectively) distinguish the mature forms of 
M. truncatula plastid-located GS from the cytosolic counter- 
parts. The dicotyledonous enzymes are also strikingly similar 
to the only structurally characterized monocotyledonous 
homologue, Z. mays GSII-la (Unno et al, 2006), with which 
they share 62-88% identical residues, including strict conser- 
vation of all of the amino acids that directly interact with the 
substrates (Fig. la). 

The susceptibihty of each of the four M. truncatula GSII 
isoenzymes to inhibition by phosphinothricin (PPT) and 
methionine-5-sulfoximine (MSG) was evaluated (Fig. lb). The 
two compounds are glutamate analogues and have been 
shown to compete with this amino acid for binding to the 
active site (Ronzio et al, 1969; Colanduoni & Villafranca, 
1986; Manderscheid & Wild, 1986). PPT and MSG are phos- 
phorylated in the presence of ATP, leading to the formation of 
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a compound similar to tiie tetrahedral glutamate-adduct 
intermediate, resulting in irreversible inhibition of the enzyme 
(Liaw & Eisenberg, 1994). The percentage of GS inhibition 
caused by different concentrations of PPT or MSO was 
determined under standard GS assay conditions. The 
concentration of inhibitor reducing the GS activity by 50% 
(IC50) was estimated for PPT as 0.76, 0.67, 0.26 and 0.72 mM 
for GSII-la, GSII-lb, GSII-2a and GSII-2b, respectively. For 
MSO titration (using 20 mM glutamate, given that the four 
GS isoenzymes were poorly inhibited under the standard GS 
assay conditions), the determined IC50 was 9.5, 2.2, 3.6 and 
428.6 mM for GSII-la, GSII-lb, GSII-2a and GSII-2b, 



respectively. Clearly, the four M. truncatula isoenzymes are 
more susceptible to inhibition by PPT than by MSO. Inter- 
estingly, GSIT2b was found to be extremely resistant to MSO 
inhibition. 



3.2. The atomic structure of MtGSII-la reveals a 
homodecameric quaternary arrangement 

The three-dimensional structure of M. truncatula GSII-la 
was determined at 2.35 A resolution by X-ray crystallography 
(Table 1) using synchrotron radiation (Seabra a/., 2009). The 
crystallographic model revealed a decameric arrangement of 





Figure 2 

Three-dimensional structure of M. truncatula GSII-la at atomic resolution, {a) Top view of the decameric assembly (D5 symmetry) along the fivefold 
symmetry axis (black line). M. truncatula GSII-la subunits establishing inter-ring contacts are coloured in the same hue: green (monomers A and F), 
orange (monomers B and G and monomers D and /) or dark blue (monomers C and H and monomers E and /). The red lines indicate the positions of the 
five twofold axes, {b) Side view of the M. truncatula GSII-la model resulting from a 90° rotation around x relative to the view in {a). Subunits are colour- 
coded as in (a), and the red and black lines represent the twofold and fivefold axes, respectively, (c) Close-up view of the inter-ring contact. In the left 
panel, M. truncatula GSII-la subunits are coloured green (subunits A and F). The corresponding Z. mays GSII-la subunits (subunits A and /) in the right 
panel are coloured red. The aromatic side chains involved in both interactions are labelled. In both molecules, a thicker ribbon represents subunit A. 
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the enzyme with D5 symmetry, with two stacked (face-to-face) 
pentameric rings. The pentameric rings of MtGSII-la are 
~130 A wide and stack to form a ~100 A tall prism-like 
structure (Figs. 2a and 2b). The ten MtGSII-la monomers 
adopt very similar conformations, as highlighted by the r.m.s.d. 
of 0.35-0.65 A for 291-311 aligned C" atoms. As expected 
from the considerable sequence conservation (Fig. la), 
M. truncatula GSII-la is also structurally very similar to the 
Z. mays homologue (Unno et al, 2006), with the superposed A 
subunits of both molecules displaying an r.m.s.d. of 1.00 A for 
321 aligned C atoms. In spite of this, however, six regions 
clearly deviate in these structures (Fig. la, pink lines). Three 



of these structurally divergent sections are located at the inter- 
ring interface (Aspl38-Trpl41, Trpl45-Tyrl57 and Pro236- 
Ala246; M. truncatula numbering used throughout), while two 
others stretch into the active site (Ala63-Ser68 and Asn302- 
Val315). Also in the neighbourhood of the active site, the highly 
conserved segment encompassing Gly285-Ile301 is disordered 
in the crystallographic model of M. truncatula GSII-la. 

The extensive contacts established between the subunits 
that compose each ring result in a buried inter-monomer 
surface of approximately 3500 A^. Additionally, each MtGSII- 
la monomer establishes fairly limited interactions with just a 
single partner from the opposite ring. The overall buried 




Nano-electrospray mass spectra of M. truncatula GSII-2a. (a) MS spectrum acquired for the intact MtGSII-2a assembly with a measured mass of 441 814 
± 75 Da, which is in good agreement with it existing as a decamer composed of ten copies of the M. truncatula GSII-2a monomer {b). (c) The spectrum 
recorded under harsh MS conditions is dominated by two charge-state series corresponding to the M. truncatula GSII-2a decamer (series centred at 41+) 

and pentamer (series centred at 27+), respectively, suggesting that the intact complex consists of two labile interacting pentamers. {d) Expansion of the 
miz region below an miz of 2500 in (c) showing M. truncatula GSII-2a monomers dissociated both from pentameric and decameric complexes in the gas 
phase. 
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inter-ring surface amounts to ~2500 A^, i.e. approximately 
500 per individual monomer-monomer interaction, high- 
lighting the weakness of the inter-ring contacts. Indeed, these 
contacts involve only the Trpl41-Tyrl50 segment of each of 
the intervening molecules (Fig. 2c). Considering the 
contacting monomers A and F from opposed rings, hydro- 
phobic interactions exist between the side chain of TrpMl'* 
and both Ilel47^ and Glyl48^, as well as between Glyl48^ and 
the side chain of Trpl41^. Polar contacts are established 
between Trpl45^ NEl and Ilel47^ O and between 0^148^* O 
and Trpl45^NEl. Finally, an additional hydrogen bond is 
formed between TyrlSO^ N and Glyl49''0 (Fig. 2c). As a 
result of this particular arrangement of the inter-ring contacts, 
a 45° rotation around the fivefold axis of the molecules 
separates the chains of any given pair of contacting monomers 
(Figs. 2a and lb). 

The active sites of MtGSII-la are located close to the 
intermonomer contacts within each oligomeric ring, similar 
to other glutamine synthetases. In this particular case, the 
structure of the apo form of the enzyme was determined and 
therefore there are no ligands present in any of the ten active 
centres. The segments Trpl25-Glul31, Serl87-Glul92, 
Glul99-Pro204, Gly245-Ser253 and Tyr328-Arg332 delimit 
the funnel-shaped active-site cavity. The strictly conserved 
amino acids Glul29, Glul31, Glul92, Glul99, His249 and 
Glu330 (Fig. la) are positioned approximately in the middle of 
this cavity, favourably located to mediate cation coordination. 
In all of the MtGSII-la monomers a segment comprising the 
Thr293-Ala299 region, homologous to the Glu-binding loop 
of bacterial glutamine synthetases (Liaw & Eisenberg, 1994; 
Wray & Fisher, 2010), is disordered, leaving open the entrance 
to the glutamate-binding site at the narrower end of the 
active-site cavity. 

3.3. MtGSII-2a is a homodecamer composed of two 
superposed pentameric rings 

The degree of sequence similarity of both MtGSIT2a and 
MtGSII-2b to other homodecameric GSIIs (Fig. la) suggests 



that these plastid-located enzymes are also decameric. In 
order to confirm the oligomeric state of MtGSII-2a in solution, 
we recorded mass spectra of this protein under nondenaturing 
conditions. A single charge-state series at an miz of around 
10 000 (Fig. 3fl) dominates the spectrum with a measured mass 
of 441 814 ± 75 Da, in good agreement with the theoretical 
value of 441 235 Da for a decameric arrangement calculated 
from the monomer mass measured by LC-MS (Fig. 35). The 
slight increase in measured mass arises owing to incomplete 
desolvation of the assembly, a common phenomenon attrib- 
uted to the retention of water and buffer molecules as a result 
of the 'soft' ionization methods required to maintain the 
integrity of the complex (McKay et al, 2006). Interestingly, MS 
spectra acquired under high activation conditions revealed an 
additional charge-state series at a lower mIz region to the 
intact complex (Fig. 3c). The mass measured for this series is 
221 925 ± 10 Da, consistent with the presence of isolated 
pentamers under the stringent conditions under which mass 
spectra were recorded. Since all decameric GSII enzymes are 
suggested to organize as two pentameric rings (Unno et al, 
2006), this dissociation behaviour of MtGSII-2a in the gas 
phase has implications for the nature and strength of its 
subunit interactions, with stable bonding within each penta- 
meric ring and a more labile and dynamic association between 
the double stacked rings (Benesch et al, 2007; Hernandez & 
Robinson, 2007). The less intense series below an mIz of 2500 
(Figs. 3c and ?>d) has a measured mass of 44 032 ± 51 Da, 
corresponding to monomeric MtGSII-2a dissociated from 
both MtGSII-2a pentamers and decamers in the gas phase. 
Overall, MS analysis of MtGSII-2a strongly supports a 
homodecameric organization of this enzyme in solution built 
as two stacked copies of a homopentameric ring. 

3.4. Swinging motion of MtGSII-2a rings 

Although the crystal structures of several GSII-la enzymes 
have been obtained, all attempts to crystallize MtGSII-2a 
have so far been unsuccessful. To gain insight into the three- 
dimensional structure of this enzyme, we collected cryo-EM 




Raw Reference-free two-dimensional averages Scheme 




Electron cryomicroscopy images of M. truncatula GSII-2a. (a) EM field showing a homogeneous distribution of particles on the ice. {b) Comparison of 
the three major classes of particles observed in the data set, showing examples of single raw images (left) and reference-free two-dimensional class 
averages (middle). A schematic representation of the averages is shown for clarity (right). 
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id) 

Figure 5 

Electron cryomicroscopy reconstruction of M. truncatula GSII-2a. (a) FSC showing a resolution of 20 A. (b) Pairs of class averages (upper row) and 
projections of the final three-dimensional reconstruction (bottom row), (c) Structure of M. truncatula GSII-2a, with approximate dimensions and 
positions of the symmetry axes (arrowed lines and pentagon for twofold and fivefold axes, respectively) defining the D5 symmetry. The numbers 1-5 and 
the numbers in parentheses indicate the positions of the five high-density blobs in the top and bottom rings, respectively. A scheme shows the relative 
rotation of the pentagons defined by the rings, (rf) Fitting of the crystal structure of M. truncatula GSII-la reported here, in which each monomer is 
coloured differently. 
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images under quasi-native conditions (Fig. 4a). More than 
7000 individual particles were selected and classified using 
reference-free two-dimensional alignment algorithms as 
implemented in EM AN (Ludtke et al, 1999) and Xmipp 
(Scheres et ah, 2005). Two major classes were readily apparent 
in the data set, corresponding to the top and side views of the 
molecule, which are related by a 90° rotation (Fig. 45). In the 
top view, the pentameric ring expected from our native mass 
spectrometry results can be observed, although the symmetry 
is somewhat difficult to appreciate from two-dimensional 
averages (Fig. Ab, top row). This is probably owing to the fact 
that the two rings are not perfectly aligned with respect to 
each other, as observed for all other GSII enzymes with 
known atomic structure. In the side view, two major and 
parallel masses of density can be identified which are 
connected by a central thinner mass corresponding to a weak 
inter-ring region (Fig. Ab, central row). These views correlate 
well with the published structures of other GSII enzymes 
(Unno et al, 2006; Krajewski et al, 2008; He et al, 2009) and 
also with the crystal structure of MtGSII-la reported here 
(Figs. 2a and 2b). 

Interestingly, we observed additional class averages in 
which the two masses of density that we attribute to the 
pentameric rings are no longer parallel (Fig. 4c, bottom row). 
The sizes of the rings in these classes are identical to those 
observed in the standard side views, but the rings come as 
close as 20 A on one edge, whereas they separate by up to 
50 A at the opposite end, generating a maximum tilt of 20°. In 
these averages, the central and weak mass observed in cano- 
nical side views is more prominent and is no longer located 
in the perfect centre of the ring mass, but rather condenses 
towards the edge where the two rings are more distant, 
allowing an approximation of the rings on the opposite edge. 
This swinging motion of the rings with respect to each other 
is likely to be owing to flexibility of the inter-ring connections. 
The presence of distinct subclasses in which the two rings are 
perfectly defined (Fig. Ab, bottom row) suggests that non- 
parallel ring particles may reflect intermediate assembly stages 
in the formation/disruption of the enzyme decamer. 



3.5. Electron cryomicroscopy structure of MtGSII-2 

To deepen our understanding of the MtGSII-2a structure, 
we applied single -particle reconstruction methods to generate 
a three-dimensional model using the collected cryo-EM 
images. Because the mass spectra of MtGSII-2a strongly 
support a double pentameric organization, fivefold symmetry 
was applied during reconstruction procedures, using only 
particles in class averages corresponding to canonical side 
views {i.e. with parallel rings). Thus, the data set from which 
the reconstruction was calculated contains 3117 particles with 
fivefold symmetry. The final volume has a resolution of 20 A 
(Fig. 5a) based on Fourier shell correlation with a cutoff value 
of 0.5 and matches the experimental data well (Fig. 5b). It 
represents the first structural information on a plastid-located 
GS enzyme. 



The structure of MtGSII-2a has an overall cylindrical shape, 
with approximate dimensions of 100 A in height and 130 A in 
diameter (Fig. 5b). It is formed by two parallel and symmetric 
rings, each composed of five monomers of the enzyme, as 
suggested by the five masses of globular density composing 
them (numbered 1-5 in Fig. 5c). When observed from the top 
{i.e. along the fivefold symmetry axis), the pentagons defined 
by the two rings are rotated by ~45° with respect to each 
other, which may have functional implications. Each of the 
rings is about 40 A thick, leaving an inter-ring space of about 
35 A that is crossed at the centre by five threads of density 
connecting the rings. These threads are only 12 A thick, which 
suggests that they could be flexible and thus rearrangements 
within them may provoke movement of the rings relative to 
each other. Fitting of the crystal structure of MtGSII-la shows 
that the overall structure is very similar (Fig. 5d), as expected 
from sequence alignment (Fig. ia). Similar results were 
obtained when the ZmGSII-la crystal structure was fitted (not 
shown). Two distinct regions, comprising about 15% of the 
total enzyme mass, appear to protrude out of the density at the 
inter-ring space, likely owing to flexibility. The first region, 
comprising three loops (Ala63-Gly68, Pro236-Gly247 and 
Gly285-Asn302), is involved in interaction with the substrates, 
and thus these loops may only become ordered in their 
presence. The second region, comprising two loops (Glnl36- 
Trpl41 and Trpl45-Glyl55), forms the inter-ring connections 
that in our reconstruction seem to adopt a shghtly different 
orientation to that observed in the crystal structure. Interest- 
ingly, both regions correspond to the areas that present 
different conformations in the two available crystal structures 
of plant GSII-la. Although the resolution of our reconstruc- 
tion does not allow atomic modelling of the active site, strict 
conservation of the residues involved in substrate recognition 
suggests that the structure of this region wiU be very similar to 
that observed for ZmGSII-la, as discussed for GSII-la. 

4. Discussion 

In this report, we present the structures of two isoforms of 
M. truncatula GSII from different subcellular compartments 
determined using complementary experimental techniques. 
Cytosolic MtGSII-la was characterized by X-ray crystaUo- 
graphy and represents the first atomic structure of a dicoty- 
ledonous plant GSII enzyme, whereas the three-dimensional 
reconstruction of MtGSII-2a obtained by electron cryo- 
microscopy provides the first structural information on a 
plastid-located GSII. The two enzymes share a common 
decameric architecture, similar to the quaternary arrangement 
of other eukaryotic glutamine synthetases (Unno et al, 2006; 
Krajewski et al, 2008; He et al, 2009), with two stacked anti- 
parallel pentameric rings that interact through a limited 
surface region forming thin inter-ring connectors. 

Despite the remarkable structural similarity between 
MtGSII-la (r.m.s.d. of 0.98-1.03 A for 320-322 aligned C" 
atoms) and its closest structurally characterized homologue, 
ZmGSII-la (Unno et al, 2006), there are some critical 
differences between the two enzymes. Since MtGSII-la was 
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crystallized in the absence of substrates or inhibitors, a 
somewhat large region (Gly285-Ile301) that comprises the 
Glu-binding loop (Thr293-Ala299) is disordered. There is also 
a concerted movement of several segments of MtGSITla, 
notably the ~180° rotation of Ala246, with its side chain now 
protruding into the glutamate-binding pocket, and the ~90° 
rotation of the upstream Trp243 side chain that packs against 
Pro248 and leads to a rearrangement of the Lys237-Asp242 
segment. TyrlSO now occupies the space left vacant by the 
Trp243 side chain, with concomitant reorganization of the 
inter-ring contact. Despite the near-perfect sequence conser- 
vation of the Trpl41-Glyl52 segment involved in establishing 
the interactions between the two pentameric rings of both 
enzymes (Fig. la), the arrangement of this surface loop is 
dramatically different from ZmGSII-la (Fig. 2c). As a 
consequence, there is a difference of ~15° in the relative 
rotation of the two superposed rings between the free 
MtGSII-la and the inhibited ZmGSII-la structures. It is 
therefore tempting to speculate that substrate binding and 
release by eukaryotic GS is connected to the rotation of the 
two homopentameric rings of the enzyme through the 
concerted movement of a limited set of regions in the vicinity 
of the active-site cleft and the inter-ring contact, and is in good 
agreement with the observed positive cooperativity of all 
isoforms of MtGSII (Seabra et al, 2013). 

The plasticity of the inter-ring contact is further under- 
scored by a previously unidentified movement of the penta- 
meric rings with respect to each other, as observed in cryo-EM 
images of MtGSIT2a. In this swinging motion, the rings almost 
contact on one side, breaking the perfect parallelism observed 
in all of the described crystal structures. We postulate that this 




Figure 6 

A small number of discrete amino-acid replacements give M. truncatula 
GSn-2b significant resistance to herbicide inhibition. The amino-acid 
sequence positions unique to M. truncatula GSII-2b and located in the 
vicinity of the MSO binding site are highlighted (stick representation) in 
the experimental model of M. truncatula GSII-la. The active site at the 
interface between subunits A and B (colour code as in Fig. 2) is depicted 
and a methionine-5-sulfoximine molecule is shown in space-filling 
representation, as found in the complex with Z. mays GSII-la (PDB 
entry 2d3b; Unno et al, 2006). 



phenomenon may be important for decamer assembly from 
two pentameric rings, as several different relative divergences 
of the rings are observed (Fig. 45, bottom row). This is possible 
because of rearrangements occurring at the level of the rela- 
tively weak inter-ring contacts, which are used as a hinge, as 
suggested above to explain the difference between the GSII- 
la structures of M. truncatula and Z. mays. The fragility of the 
inter-ring contacts is supported by native mass-spectrometric 
data, showing that stringent mild conditions provoke the 
dissociation of decamers into pentamers. Recently, active 
pentameric GS rings have been detected by native gel elec- 
trophoresis in root nodules of M. truncatula (Seabra et al, 
2013), and GS dissociated active rings have also been reported 
to occur in the roots of Beta vulgaris (Mack, 1998). According 
to the structural features reported here, there is no structural 
impairment of the preservation of catalytic activity in GS 
pentamers. However, differences in catalytic activity between 
the two oligomeric forms are predictable and could represent 
an additional mechanism of GS activity regulation in plant 
cells. The finding that active GS pentamers are exclusively 
detected in root nodules suggests differences in the assembly 
and/or stability of the enzyme in this particular organ, which is 
of particular interest. 

The availability of the three-dimensional structures of 
closely related GS isoenzymes may be a valuable tool for the 
design of new herbicides. The four MtGSII isoenzymes differ 
in their catalytic properties (Seabra et al, 2013) and, as shown 
here, in their susceptibility to inhibition by PPT or MSO. The 
IC50 values that we have determined for both inhibitors are in 
accordance with previous studies showing that, overall, plant 
GSII is 5-10 times more effectively inhibited by PPT than by 
MSO (Evstigneeva et al, 2003). Nevertheless, while the four 
MtGSII isoenzymes showed a similar susceptibihty to inhibi- 
tion by PPT, their sensitivity to MSO is strikingly different 
(Fig. lb). The surprisingly elevated resistance of the MtGSIT 
2b isoenzyme to MSO (ICso of 428.6 mM) when compared 
with GSII-la, GSII-lb and GSII-2a (IC50S of 9.5, 2.2 and 
3.6 mM, respectively) is noteworthy. 

This remarkable functional difference is likely to arise from 
unique amino-acid changes with respect to the other three 
isoforms, which in M. truncatula GSIT2b are limited to a few 
positions (Fig. la). The majority of these substitutions are 
located on the enzyme surface and at a considerable distance 
from the active site, most likely being rather neutral. However, 
a few unique amino-acid changes occur in the vicinity of the 
MSO binding site and are therefore hkely to modulate the 
sensitivity of M. truncatula GSIT2b to herbicides (Fig. 6). 
Notably, the Ile364-Gly365-Lys366 triplet (Val-Ala-Asn in all 
other M. truncatula GSII isoforms) precedes the invariable 
Arg367 (Arg311) involved in substrate binding, which in turn 
could estabhsh contacts with the neighbouring Thrl24 (Ser68). 
In close proximity, the strictly conserved Arg372 (Arg316), 
also a substrate ligand, leans against the Ile364-Gly365-Lys366 
segment. Thus, its stabilization might be affected by the 
increased flexibility brought about by the replacement of an 
alanine by a glycine at position 365. Also, the metal-coordi- 
nating Glu248 (Glul92) is preceded by an alanine residue in 
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M. truncatula GSII-2b, instead of the more flexible glycine 
found at this position in all other isoforms. Finally, and most 
strikingly, at one of the hinges of the glutamate-binding loop 
Asp358 (Asn302) is in close proximity to the invariable His333 
(His277) and Lys334 (Lys278) side chains, and could therefore 
be involved in stabilizing the open conformation of the 
glutamate flap by estabUshing electrostatic interactions with 
these residues. The flexibility of the glutamate flap is essential 
for the catalytic activity of the enzyme, allowing the exchange 
of newly synthesized glutamine for glutamate. In its closed 
conformation, the glutamate-binding loop shields the reaction 
intermediates from solvent and harbours the essential Glu353 
(Glu297) that binds the transition-state analogue arising from 
MSO phosphorylation (Liaw & Eisenberg, 1994). Therefore, 
stabilization of the open conformation of the glutamate- 
binding loop is expected to result in increased resistance to 
inhibition by MSO, as was observed in mutational studies of 
Bacillus subtilis OS, where changes in the glutamate flap led to 
up to 200-fold increased levels of resistance to inhibition by 
MSO (Wray & Fisher, 2010). 

In conclusion, we have unveiled the quaternary architecture 
of dicotyledonous glutamine synthetase, and report a common 
decameric arrangement of cytoplasmic and plastid-located 
enzymes from the model legume M. truncatula. Both enzymes 
are composed of two juxtaposed homopentameric rings linked 
by a surprisingly flexible inter-ring contact. Moreover, we 
propose a zipper-like mechanism for ring-ring assembly, 
dependent on rearrangements of the inter-ring loops. Despite 
the structural similarity of the reported enzymes, there are 
striking differences in herbicide resistance among the distinct 
M. truncatula GS isoforms. The structural models now 
reported allowed the identification of a limited set of discrete 
and unique amino-acid substitutions in MtGSll-2b that 
account for its distinctive MSO resistance and pave the way 
for the design of herbicide-resistant plants. 
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